OPEN 3 ACCESS Freely available online 



•0-PLOS I ONE 



Determinants of Genetic Structure in a Nonequilibrium tm\ 
Metapopulation of the Plant Silene latifolia qzssp- 

Peter D. Fields*", Douglas R. Taylor 

Department of Biology, University of Virginia, Charlottesville, Virginia, United States of America 

Abstract 

Population genetic differentiation will be influenced by the demographic history of populations, opportunities for migration 
among neighboring demes and founder effects associated with repeated extinction and recolonization. In natural 
populations, these factors are expected to interact with each other and their magnitudes will vary depending on the spatial 
distribution and age structure of local demes. Although each of these effects has been individually identified as important in 
structuring genetic variance, their relative magnitude is seldom estimated in nature. We conducted a population genetic 
analysis in a metapopulation of the angiosperm, Silene latifolia, from which we had more than 20 years of data on the spatial 
distribution, demographic history, and extinction and colonization of demes. We used hierarchical Bayesian methods to 
disentangle which features of the populations contributed to among population variation in allele frequencies, including 
the magnitude and direction of their effects. We show that population age, long-term size and degree of connectivity all 
combine to affect the distribution of genetic variance; small, recently-founded, isolated populations contributed most to 
increase F S j in the metapopulation. However, the effects of population size and population age are best understood as 
being modulated through the effects of connectivity to other extant populations, i.e. F ST diminishes as populations age, but 
at a rate that depends how isolated the population is. These spatial and temporal correlates of population structure give 
insight into how migration, founder effect and within-deme genetic drift have combined to enhance and restrict genetic 
divergence in a natural metapopulation. 
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Introduction 

Sewall Wright [1] introduced the study of geographic popula- 
tion structure and its genetic consequences with the Island Model. 
In the Island Model, populations are assumed to be stable and 
interconnected by migration, with population differentiation (jFst) 
at neutral loci generated by drift and diminished by gene flow 
among populations. When the simplifying assumptions of the 
island model are violated, a number of additional factors are 
necessary to understand population genetic differentiation [2,3] . In 
metapopulations, where demes experience extinction and recolo- 
nization, there are broad conditions under which the founder 
effect can be a powerful structuring mechanism [4], and 
population differentiation can occur despite high levels of gene 
flow. The importance of founder effects depends on several 
parameters including migration, extinction and colonization rates, 
the number of founding propagules, as well as the fraction of 
demes from which the colonists originate [2,3]. 

Empirical evidence that founder effects can drive population 
differentiation has been gathered from age-structured populations. 
For example, elevated population structure among newly estab- 
lished demes in metapopulations is consistent with founder effects 
generating increased population structure, with subsequent gene 
flow reducing genetic differentiation among populations as they 
age [5-9]. Although these classic studies demonstrate that 



metapopulation processes (extinction and recolonization) can 
influence genetic differentiation, it is important to estimate their 
importance relative to other structuring mechanisms such as local 
adaptation, which will also contributed to increased population 
genetic structure [10,11]. 

Recent advances in evaluating the causes and consequences of 
population structure have benefitted from the development of a 
Bayesian model-based approach known as the F-model [12-17]. 
This approach has an advantage over previous methods in that 
rather than estimating "global" values of genetic differentiation, it 
takes into account the biological reality that local populations 
differ in their effective sizes and migration rates and, therefore, 
their degree of differentiation from other populations [12]. Foil 
and Gaggiotti [ 1 8] introduced a hierarchical formulation of the F- 
model that uses population-specific measurements to obtain priors 
for F &T , then estimates the proportional contribution of popula- 
tion-specific variables generating among population variance in 
allele frequencies. 

We used the hierarchical Bayesian approach of Foil and 
Gaggiotti [18] to understand variation the causes Fgx m a well 
studied plant metapopulation. Our objective was to estimate how 
population structure results from equilibrium process such as those 
in the island model (e.g. drift, migration) versus non-equilibrium 
factors that operate in metapopulations (e.g. extinction and 
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recolonization with founder effects). To accomplish this, we 
applied the F-model formulation to a long-term data set of the 
plant, Silene lalifolia, where the size and spatial distribution of 
demes was known, and could be combined with information about 
their demographic history and age structure. We used fine-scale 
sampling, microsatellite genotyping, and the application of an F- 
model based approach to explore the spatial and temporal features 
of populations that affect their contributions to the observed 
distribution of F^t, including simultaneous estimates of the 
magnitude and direction of these effects. 

Materials and Methods 

Study Organism 

Silene lalifolia Poir. (= S. alba, Caryophyllaceae) is a short- 
lived perennial plant that is broadly used as a model system for 
studying sex determination, sex chromosome evolution, host- 
pathogen dynamics, biological invasions, organelle evolution, 
sexual dimorphism, sex ratio evolution, and evolution in 
structured populations [19]. 

We studied a S. lalifolia metapopulation located in Giles and 
Craig Counties, Virginia, USA (Figure 1). This region has been 
the subject of a 20+ year study of population dynamics and genetic 
structure in more than 800 populations [5,20-25]. An annual 
census of an approximately 25 x25-km area adjacent to Mountain 
Lake Biological Station has been conducted since 1988, recording 
the location and sex of S. lalifolia individuals along ~ 150 km of 
predominantly roadside habitat. The structure of the data and 
how it was collected is reported in Antonovics et al. [20] . Briefly, 
the roadside habitat is divided up into ~40 m grid units which we 
classify as populations. Previous studies have indicated that the 
utilized grid system may encompass one or two genetic neighbor- 
hoods on average [26,27], levels of extinction and colonization 
scale consistently with grid segment size [20,28], and that dispersal 
distances are such that populations separated by 20 meters were 
almost panmictic, while populations separated by 80 meters were 
nearly isolated [24]. Dispersal of pollen is insect mediated, 
primarily noctuid moths, and seed dispersal is passive and happens 
at a much smaller scale [22]. Much longer distance dispersal is 
likely mediated by human activity [29], with road grading and 
mowing as possible mechanisms [24]. Which units are occupied 
and the number of plants in each occupied population are 
recorded annually. Important phase transitions, such as extinction 
and colonization events, are confirmed with a second census 
during the same season. Time since colonization (or population 
age) is based on the year plants were first observed in a given grid 
unit. We identified extinction as the disappearance of plants from 
a grid unit for a single year. A previous study has indicated that 
seeds can remain viable in the soil for approximately four years, 
raising the possibility of a seed bank [30]. For our purposes, we 
treat the recolonization of a population by propagules in the seed 
bank as roughly equivalent as recolonization by propagules from a 
nearby site. The census data provide the demographic data, spatial 
relationships among populations, and extinction/recolonization 
dynamics used in this study (Table 1). 

Spatiotemporal Data 

We estimated characteristics of populations that are thought to 
significandy influence the magnitude of genetic differentiation 
among demes in a metapopulation, including population size, 
population connectivity (a characteristic that will influence gene 
flow between populations), and population age (a factor that 
estimates the when founder effects occurred and opportunities for 
subsequent migration). 




Figure 1. Map of the focal populations of the S. latifolia 
metapopulation located in Giles and Craig County, VA, USA 
sampled in the presented analysis. Circles represent individual 
populations, where the size of the circle indicates the total number of 
plants located within our grid. Black lines represent both the assumed 
grid of population arrangement and small country roads. Due to the 
topology of the focal area (mountain, valley systems), it was assumed 
that pollinators move along the linear grid, rather than crossing over 
ridges. 

doi:1 0.1 371/joumal.pone.01 04575.g001 

Given our incomplete knowledge of the contribution and life- 
span of an individual population's seed bank, we estimated 
harmonic mean population size (hereafter population size) based 
upon the number of plants (both flowering and vegetative) 
occupying a given grid-unit throughout the time of the annual 
census as: 

Population size estimates derive from the occupancy of a given 
site throughout the full census period. Detailed surveys on a subset 
of the metapopulation have been conducted to estimate the extent 
that our non-invasive census protocol sacrifices precision and 
accuracy. In the survey, we carefully searched the vegetation and 
when plants were found, the shoots were carefully traced down to 
the ground to distinguish individuals that were rooted close 
together. Non-flowering plants were also counted. Population sizes 
in the detailed survey were highly correlated with population size 
from the metapopulation census as assessed with a Pearson 
Product Moment Correlation (0.74, P<0.001). Population size in 
the metapopulation census (ft = 14.32) is lower than in the detailed 
survey (p =36.91) because clustered plants many be counted as 
single individuals in the annual census. The census, therefore, 
accurately measures relative population sizes, but likely underes- 
timates absolute size. 
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To estimate population age, we assumed sites unoccupied for a 
single year went extinct, and that plants occupying that site in the 
future were classified as colonizations. This reflects the conserva- 
tive assumption that recolonization from neighboring sites and 
recolonization from the seed bank involve similar bottlenecks of 
genetic diversity; if this assumption were violated, and new 
colonization from a seed bank did not involve the same degree of 
bottleneck, this would tend to obscure any differences we observed 
between younger and older populations. Following extinction and 
re-colonization, age is calculated as the number of years a site has 
been occupied up until the time of collection. Given this 
operational definition of population age, individual sites could 
vary in age from one to 2 1 years (the extent of the demographic 
census started in 1988). Populations in the age class of 21 years 
might reasonably be considered as a heterogeneous grouping of 
extant populations given that some populations may actually be 
older than the extent of the current long-term census. 

There are a number of methods described in the literature that 
have been used to estimate population connectivity. The majority 
of ecological studies have used a nearest neighbor/patch 
approach, or distance to multiple neighbors within a limited 
neighborhood of a focal patch (or buffer) [31]. However, these 
measures have been shown to be poor predictors of ecologically 
important metapopulation dynamics such as colonization poten- 
tial. To measure population connectivity, we generated two 
composite variables, each of which includes a negative exponential 
dispersal kernel and accounts for distance to all other potential 
extant populations for gene flow [32]. Each variable has the 
following structure: 

C,-= ^/> 7 -exp(-a<%),(l) 

where pj is in C^ ge the proportion of the censuses for which a 
given site was occupied (1/21 - 1) or in Cy, the population size 
over the life span of the population, of some population j; a is a 
parameter scaling the effect of distance on dispersal, and dy is the 
road network based distance between a target population i and 
source population j [31,32]. Pair-wise distances were calculated 
using a network constructed based upon the public roadway 
system, using ArcGIS (ESRI) Network Analyst tool. Given the 
mountain-valley geographic topology of the area, this network- 
based approach is more appropriate than standard Euclidean 
distances. C^ g( , or Cjf are indicative of a population's probability 
of receiving migrants, whether through seeds or pollen, though 
these are modulated by the opportunity for gene flow over-time or 
through apparency to pollinators, respectively. 

Sampling 

We focused on one of the nine metapopulation sections, which 
are classified with an ad hoc partitioning of the mountain-valley 
system topology of the area, and is consistent with the expectation 
that each are mostly genetically isolated [27]. We sampled plants 
from 33 spatially distinct populations during peak flowering in the 
summer of 2008. Sites occupied by < 4 plants were excluded since 
we could not reliably estimate gene frequencies in our estimation 
of F ST . S. latifolia is identified as in introduced weed in the study 
area, and does not require specific permits to collect tissue 
samples. We collected leaf tissue from every plant in the 
population, or up to 50 individuals in the larger populations, 
and stored the leaves with silica gel (Sigma, USA). 

Genomic DNA was extracted and amplified following estab- 
lished microsatellite techniques for S. latifolia [33]. We geno typed 



each plant at 8 microsatellite loci. Microsatellites were derived 
from multiple sources [34-36]. PCR amplification was conducted 
using previously published methods [37]. Briefly, PCR products 
were amplified with the forward primer end-labeled with a 
fluorescent dye, either 5(or 6)-FAM, NED, TAMRA, JOE, or 
VIG Three to four PCR products of different loci were then 
pooled together and added to a loading buffer containing 
formamide and GENESCAN 400HD ROX size standard 
(Applied Biosystems, USA). Following five minutes of denaturing 
at 95°C, fluorescendy labeled fragments were separated on an 
Applied Biosystems 3130 sequencer and analyzed with GENE- 
MAPPER v3.0 software (Applied Biosystems, USA). Allele binning 
was accomplished using the software TANDEM [38] . 

Statistical Methods 

Overall, the data involved 730 plants each with a multi-locus 
genotype derived from eight microsatellite markers. The plants 
were associated with one of 33 populations that had size, age and 
connectivity data collected from the long-term demographic 
dataset. 

We used the program GenoDive [39] to estimate global summary 
statistics of population structure for our molecular markers, 
including observed (Ho) and expected (H E ) heterozygosity, F ST 
andjost's D [40] . We used the the analytical software R [4 1] and the 
package hierfstat v. 0.04-10 [42] to calculate rarefied allelic richness. 
We used the hierarchical Bayesian method of Foil and Gaggiotti 
[18], implemented in the program GESTE v. 2, to evaluate the 
effect of spatial and temporal characteristics of populations on the 
magnitude of genetic differentiation among populations. The 
method implemented in GESTE estimates -Fst values for each 
sampled population using the approach first proposed by [14] and 
relates them to environmental factors using a generalized linear 

model framework [43] . Specifically, In (^st j \— Fg T ) ' s related to 

models of the form ao+ / J ._ 1 a iGji, where ctois a constant, a,is a 
regression coefficient of the effect of factor Gyfrom some pop j 
affecting the response variable, in this case individual population 
F S t values. We considered four factors, thus generating 16 
alternative regression models The simplest null model included 
only a constant term, whereas the most complex model included a 
constant and all four factors. Because the two connectivity variables, 
CU ?c and Cjj, were found to be significantly correlated with a 
Pearson product-moment correlation coefficient of 0.6141 (95% 
confidence interval (CI) 0.3432-0.797; P = 0.0001), we additionally 
ran GESTE excluding one or the other connectivity measure. 
GESTE does not allow for testing a full model that would include all 
four pairwise interactions. The method provides posterior proba- 
bilities for the alternative models using a RJ-MCMC approach [18]. 
Under this framework, the model with the highest posterior 
probability is the one that best explains the data [18,43]. We 
followed the method of Gaggiotti et al. [43], using 10 pilot runs of 
1 000 iterations to obtain parameters of the proposal distributions 
used by the MCMC, followed by an additional burn-in of 5xl() 6 
iterations and a thinning interval of 50, with a final iteration sample 
size of 60,000 on which the model fit probability was based. 
Convergence was assessed visually using the Plot GESTE.exe 
function included in the GESTE v.2 distribution. Using this 
method, we identified the model that best explained the observed 
genetic structuring. The magnitude and direction of a spatiotem- 
poral character on genetic structuring was inferred from estimates of 
the regression coefficients from the model with the highest posterior 
probability. 
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Table 2. Global estimates of genetic diversity and variation in allele frequencies. 





Locus 


N 


H„ 




D 


Fsr 




10 


0.338 


0.592 


0.158 


0.09 


5/af_J2 n 


9 


0.303 


0.673 


0.337 


0.135 


5/af_33 n 


3 


0.07 


0.13 


0.051 


0.231 


slat_48* 


2 


0.078 


0.202 


0.047 


0.143 


slat_72^ 


12 


0.325 


0.593 


0.103 


0.055 


slat_85^ 


16 


0.324 


0.538 


0.178 


0.125 




43 


0.693 


0.819 


0.245 


0.051 


SV_11 3 


11 


0.305 


0.573 


0.19 


0.122 


Overall 


13.25 


0.304 


0.515 


0.131 


0.103 



Variables are the number of alleles (A/), observed heterozygosity {H 0 ), Expected heterozygosity (H s ), Jost's D (D) and global variation in allele frequencies {F ST ). Overall, 
populations showed a high degree of substructure, as has been observed in other S/7ene metapopulations. Population structure was high for each marker. Our lowest 
F ST corresponded to the only marker composed of a dinucleotide repeat. 1 [34], 2 [35], 3 [36]. 
doi:1 0.1 371 /journal.pone.01 04575.T.002 



Results 

Global summary statistics (Table 2) revealed a high degree of 
population structure, with an overall i r sT = 0.103, and a range 
among markers of 0.055 to 0.231. This was similar to previous 
estimates in this metapopulation using allozymes [0.134, 21]. The 
observed levels of population genetic structure are also consistent 
with other published datasets on plant metapopulations using 
microsatellite markers, including other Silene species. Previous 
research has shown that the sub-structuring of alleles is not the 
consequence of PCR artifacts such as null alleles [34—36]. We 
estimated a significantly negative Pearson product-moment 
correlation coefficient between the genetic diversity of individual 
populations and that population's contribution to F ST (r = — 
0.3601, 95% CI -0.6260, -0.0191; P = 0.0396; Table SI). D had 
a mean of 0.131, and a range of 0.047-0.337. 

The magnitude and direction of the effect of each of the 
spatiotemporal characters on genetic differentiation was inferred 
from the estimates of the regression coefficients for the most 
probable model. The four parameters affected genetic differenti- 
ation in directions that were consistent with both theoretical 
expectations concerning the role of spatial-temporal dynamics 
affecting the distribution of F ST across a metapopulation [8] , as 
well previous research on the system (Table 3, 4; Fig. S1-S3). For 
example, recent founder effects and small population size 
enhanced F sx , as was previously described theoretically by 
Whitlock [8], and has been shown empirically in a number of 
other systems across a broad range of taxa known to occur as a 
metapopulations [6,9] [5,6,9], including the focal plant metapop- 
ulation. However, the parameters of population age and 
population size were not individually components of models with 



the highest posterior probabilities. Instead, a composite parameter 
that combined population age and the degree of connectivity, 
Cjge> provided the highest posterior probability (0.752) model 
(Table 3, 4; Fig. S2), consistent previous studies by Moilanen and 
Nieminen [31] for predicting extinction and colonization process- 
es. To identify which population-level character contributed the 
most to global F ST , we summed the posterior probabilities of 
models that included individual factors. Among these, population 
age (CAge > with a posterior probability of 0.994), emerged as the 
one factor that had an overwhelming effect on a population's 
contribution to the global value of F ST (Table 3). The 14 other 
models attained <0.1 posterior model probability. 

Comparison of the GESTE derived estimates of individual 
factors contribution to variation in Fst across individual models, 
i.e. all models including a given factor, showed that population age 
was the most important determinant of population genetic 
structure, followed by population size, and that each of these 
were influenced by the degree of connectivity. GESTE analyses 
run excluding population size, Cy, resulted in the Ca^ factor 
remaining the dominant factor increasing F ST (Table S2.a). The 
reverse analysis, which excluded CA ge and included Cjj, resulted 
in a model with a low posterior probability (Table S2.b) 

Posterior estimates of the regression parameters of the model 
that only included Constant + C^were significandy negative 
(p= -0.439, mode =-0.447, 95% HPDI [-0.680; -0.200]), 
indicating that the initial increase in F^t resulting from founder 
effects is reduced over time, presumably by subsequent migration 
(Table 5; Fig. S3). Posterior estimates of the regression parameters 
of the next best model indicated a consistent effect of population 
age, C Age (p,= -0.575, mode = -0.600, 95% HPDI [-0.857; - 



Table 3. Sum of posterior probabilities of models that include a given factor. 





Factor 


Sum of the posterior probabilities 


Population Age 


0.0935 


N 


0.0541 


CAge 


0.994 




0.124 



Bold value indicates factor with highest score. 
doi:l 0.1 371 /journal.pone.01 04575.t003 
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Table 4. Posterior probabilities of all 16 models. 





Model 


Pr 


Factors included 


5 


0.75 


Casc 


13 


0.10 


c c 


6 


0.08 


C Ase , Population Age 


7 


0.04 


C Age ,N 


15 


0.01 




14 


0.01 


CA ge ->Cjj, Population Age 


1 


0.00 


Null(ConstantOnly) 


8 


0.00 


C J 4 ?e) iV,Population Age 


16 


0.00 


C^ejC^A^Population Age 


2 


0.00 


Population Ase 


9 


0.00 




3 


0.00 


N 


10 


0.00 


(^Population Age 


4 


0.00 


7V,Population Age 


11 


0.00 




12 


0.00 


Cjf,N, Population Age 



Our most probable model (bolded) included the composite variable of population age and population connectivity. The second most probable model includes 
composite variables of population age and connectivity, and population size and connectivity. 
doi:1 0.1 371 /journal.pone.01 04575.t004 



0.265]), and a positive effect of populations size, Cjj (j.i = 0.207, 
mode =0.201, 95% HPDI [-0.079; 0.472]) indicating that 
younger and smaller populations contributed proportionately 
greater to the overall magnitude of genetic differentiation, and 
that proportional effects were modulated through population 
connectivity. Regression parameters for models excluding one or 
the other connectivity factor were of similar direction and 
magnitude (Table S3a, S3b). 

Discussion 

In the present study, we simultaneously identified several 
important spatiotemporal parameters that contributed to popula- 
tion genetic differentiation. While previous research on other plant 
metapopulations has shown that founder effects contribute to 
population differentiation [5-7,21,22,44-46], none of these 
previous studies have been capable of quantifying the relative 
magnitudes of these effects. 

In the Silene latifolia metapopulation, population size, degree of 
connectivity and population age were each important determi- 
nants of population structure, but composite characters (so-called 
connectivity scores; [31]) were most important. This implicates the 
extinction and recolonization of demes, followed by subsequent 
opportunities for gene flow as a population ages, as the most 
important driving force in the genetic differentiation among 
populations in this system. This point is particularly clear when 
observing Population 1 which is part of the oldest age class but also 
maintains one of the higher F^r values due to isolation (Figure 2; 
upper right quadrant). The genetic consequences of founder effects 
and within-population genetic drift, as modulated by long-term 
population size, were influenced by the degree of population 
connectivity and hence the opportunities for gene flow. 

The S. latifolia metapopulation is neither a true Island Model, 
nor an idealized metapopulation [47]; populations are character- 
ized by frequent colonizations and extinctions, on the order of 5- 



20% per year [20,27,48,49], but populations vary in size, dispersal 
is limited, and within population dynamics are important relative 
to the time scale of the study. Colonizations are likely to be source- 
size and distance dependent, as is typical in other metapopulation 
systems [50-52], though this point has not been explicitly tested. 
Our findings support the general notion that founder effects 
during colonization can enhance genetic differentiation among 
populations [5], but we additionally show that the magnitude of 
these effects are large (Figure 2) relative to the structuring 
mechanism assumed in most models, i.e. genetic drift among 
extant demes as described by Wright's Island Model [1]. 

It is a nearly universal observation that a given species range will 
be composed of populations that are patchily distributed in space. 
There are clear similarities in the extinction/colonization and 
incidence parameters in Silene versus other model systems; e.g. the 
Glanville fritillary (Melitea cinxia) metapopulation in Finland [53]. 
Spatial patterns of colonization and extinction responded similarly 
to scaling in regional studies of sunflowers in the mid-west 
(Helianthus annum; [27])(Helianthus annuus; [28]). Metapopula- 
tion dynamics have also been documented in the fresh water 
crustacean Daphnia longispina and D. magna (Cladocera), where 
the founding of new populations by single or few individuals was 
followed by predominantly clonal reproduction. This lead to an 
initially large increase in F^t in newly founded populations 
followed by the gradual dissolution of variation in allele 
frequencies through limited gene flow [9]. While the direct 
quantitative values will obviously vary from system to system, one 
might expect that non-equilibrium dynamics that result from the 
extinction and recolonization of local demes would be most 
important in systems like these, where there is a high turnover of 
demes. On the other hand, drift can be expected to be a less 
powerful structuring mechanism when demes are larger and more 
permanent. For example, Cosentino [32] show that connectivity 
and wetland area were the most important factors driving 
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Table 5. Posterior estimates of regression parameters for the model with the highest posterior probability. 



Regression coefficient Factor Mean Mode 95% HPDI 

a 0 Constant -2.26 -2.25 [-2.50; -2.00] 

a 3 G,, -0.435 -0.445 [-0.687; -0.193] 

a 1 - 0.369 0.316 [0.1 68; 0.61 7] 

Parameter estimates are consistent with theoretical expectations (e.g. older and larger populations contribute proportionately less to the global F S t). 
doi:l 0.1 371 /joumal.pone.01 04575.T.005 



variation in allele frequencies amongst populations of the tiger 
salamander Ambystoma tigrinum. 

The observed global F ST in the present study was lower than 
has been previously reported. Specifically, McCauley [21,44] 
reported mean estimates Fst — 0.134 in seven allozyme markers. 
Previous research has suggested that the F ST of microsatellite 
markers will be smaller than for allozymes due to the higher 
heterozygosity of microsatellite markers and the inverse relation- 
ship of Fst an d marker heterozygosity [54]. In addition, the 
estimated F^i in the present study is likely an underestimate due to 
our sampling scheme, which include almost all individuals rather 
than a small sample. Specifically, Whidock [8] describes the 
sampling scheme used on a metapopulation of the forked fungus 
beetle (Bolitotherus cornutus) as encompassing almost the entirety 
of the population. In such cases Whitlock [8] suggests the use of a 
hypergeometric sampling correction to account for the sampling 
error normally assumed by a binomial sampling error. Recently 
Wood el al. [55] describe a more recent analysis of the same focal 
Bolitotherus cornutus metapopulation, though not exactly the 
same section as Whitlock [8] . Wood el al. [55] report much lower 
values of FgT values for microsatellite markers, which might be 
similarly affected by the issues described by Edelaar et al. [54] 
wherein markers showing larger numbers of alleles, e.g. microsat- 
ellites, will show a general trend towards smaller values of F$i 
compared to markers that have fewer alleles (allozymes), but also 
because a multinomial-hypergeometric correction was not includ- 
ed for the estimation of F S t- 



Theoretical explorations of metapopulation dynamics have 
shown that spatiotemporal characters should have direct effects on 
each other, and we found empirical evidence for this. For example, 
we found a significant positive pairwise correlation (Pearson 
Product Moment Correlation = 0.478, P<0.01) between 
population size and population age; older populations also tended 
to be larger. Thus, while we show that both of these effects are 
potentially important in generating genetic divergence, our ability 
to fully disentangle the effects of age and population size is limited. 
Neither population size nor age was significantly correlated with 
population connectivity (P>0.05). 

While the present approach provides insights into the biology of 
the Silene metapopulation, additional insights would be gained 
from sampling populations over time. Lamy et al. [56], 
longitudinally sampled populations of the freshwater snail, 
Drepanotrema depressissimum, to show that demographic survey- 
derived assessments of extinction were incorrect; rather than a true 
metapopulation, they detected persistent demes with genetic 
diversity driven by different population sizes and rates of 
immigration. Similarly, Robinson et al. [57] applied a novel use 
of approximate Bayesian computation (ABC; [58]) to genetic data 
from a Daphnia magna metapopulation described by Haag et al. 
[9]. They concluded that some populations act as persistent 
sources, similar to a continent-island model. 

Discernment of presence/ absence of extant populations within 
a census is highly efficient in S. lalifolia (see above). However, 
previous research has also suggested the presence of a seed bank 
[30], though allele frequencies observed in the seed bank are 
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Figure 2. The effect of population age, connectivity, and population size on F sr . Individual population F ST is represented by the size of the 
circle, where larger circles represent larger F ST (A) Simultaneous effect of population connectivity and population age effects on FST, and (B) of 
population connectivity and population size on F 5T . 
doi:1 0.1 371 /journal.pone.01 04575.g002 
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highly consistent with above ground cohorts [30]. We believe our 
ad hoc measure of population extinction is consistent with the 
theoretical expectations of the effect of founder effects on the 
distribution of variation in allele frequencies [8]. Temporally 
spaced samples as in Lamy et al. [56] could more conclusively 
confirm this particular hypothesis, as well analyses such as the one 
presented by Robinson et al. [57]. Further investigations such as 
those carried out by Lamy et al. [59] and Freville et al. [60], where 
the expectation of a seed bank/ resting form is assumed as part of 
larger model of demographic dynamics, will allow for a better 
understanding misidentification of extinctions would effect our 
inference of metapopulation dynamics. 

Thus far, we have focused on evolutionary processes that are 
driven by non-selective, drift related processes. However, previous 
studies have indicated the potential for metapopulation structure 
to have a significant outcome on selective dynamics and vice versa 
[61]. Since we derived our molecular markers from non-coding 
regions, they were assumed to be neutral and to be unlinked to 
functionally important genomic regions [34-36]. Although our 
panel of markers does not exhibit discontinuities in their 
distribution characteristic of outliers, given our limited panel size 
and the complexity of detecting outliers in hierarchically 
structured populations [62], we have only limited power to test 
this assumption. 

Natural selection could have a powerful effect on population 
structure, even when the alleles under selection are not closely 
linked to marker loci. Studies in this S. latifolia metapopulation 
have confirmed that population genetic structure, where individ- 
uals in closer proximity tend to be more related than expected 
from chance, can have a significant negative average effect on 
individual fitness through the expression of deleterious recessive 
alleles [25]. Because inbreeding depression will reduce average 
fitness and therefore population size, drift processes could be 
enhanced through selective reductions in population size. Because 
many of the young, recently colonized populations will experience 
inbreeding depression, gene flow from other populations may be 
enhanced beyond neutral expectations due to heterosis [63,64]. 
This initial enhancement in gene flow could, combined with 
stochastic dynamics associated with non-equilibrium metapopula- 
tion conditions, enhance the observed reductions in population 
specific Fst over time. Thus, our analyses are unable to 
disentangle neutral and the selective effects of spatiotemporal 
metapopulation structure in generating variation in allele 
frequencies among populations. 

The present study demonstrates that the appropriate combina- 
tion of long term ecological data and population genetic analyses 
may be a powerful tool for studying the mechanisms that generate 
population structure. Spatiotemporal characters of populations 
have long been hypothesized to have significant effects in 
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